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We present results from two-photon photoassociative spectroscopy of the least-bound vibrational 
level of the X 1 E^" state of the 88 Sr2 dimer. Measurement of the binding energy allows us to determine 
the s-wave scattering length, agg = —1.4(6) ao- For the intermediate state, we use a bound level on 
the metastable 1 So- 3 Pi potential, which provides large Franck-Condon transition factors and narrow 
one-photon photoassociative lines that are advantageous for observing quantum-optical effects such 
as Autler-Townes resonance splittings. 



I. INTRODUCTION 

Precise knowledge of interactions between ultracold 
atoms has enabled spectacular advances in the produc- 
tion and study of quantum gases The most accurate 
tool for determining those interactions is spectroscopy 
of bound molecular states, such as two-photon photoas- 
sociative spectroscopy (PAS) in which two laser fields 
couple colliding atoms to a weakly bound state of the 
ground molecular potential via a near-resonant interme- 
diate state (Fig. [1]). Two-photon PAS has been used to 
measure binding energies in Li HQ, Na 0, K [jj, Rb 
1 He , and Yb [§[ . Each of these measurements pro- 
vides accurate determination of the atomic s-wave scat- 
tering length (a) and understanding of the path towards 
quantum degeneracy and behavior of resulting quantum 
fluids. Here, we report two-photon PAS of 88 Sr and deter- 
mination of a for the ground molecular potential (X 1 !]^). 
Through mass-scaling, we also determine a for all stable- 
isotope collisional combinations. 

Alkaline-earth atoms such as Sr, and atoms with 
similar electronic configuration, differ significantly from 
alkali- metal atoms that are typically used in ultracold ex- 
periments. They have a closed-shell ground state struc- 
ture, numerous isotopes including spinless bosons, and 
metastable triplet levels that lead to novel laser-cooling 
techniques @ and interactions [H EH- They present 
many new opportunities for the study and application of 
ultracold atoms, such as optical frequency standards [12j |. 
long-coherence-time interferometers 13j, and Bose and 
Fermi quantum degenerate gases and mixtures [lil . ITol ] . 
While the collisional properties of these atoms have been 
the subject of intense study [U, EE El El HI HH, H3, HI , 
until now precise scattering length values have only been 
published for Yb ||. 

For two-photon PAS of Sr, we utilize an intermedi- 
ate state that is bound in the U potential that cor- 
responds to the 1 So + 3 Pi asymptote at long range. 
The spin- forbidden So- Pi intercombination transition 
at A = 689 nm is weakly allowed due to spin-orbit cou- 
pling of the 3 Pi state with the lowest-lyin g 1 P i level 
[24j . PAS involving an intercombination line [l6l . I20L [2oT | 
transition differs qualitatively from spectroscopy near an 
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FIG. 1: Two-photon PAS diagram. The energy of two well- 
separated 1 So atoms at rest is taken as zero. E g is the kinetic 
energy of the colliding atom pair. Eb 1 is the unperturbed 
energy of the bound state of the excited molecular potential 
that is near resonance with the free-bound laser. Eb 2 (< 0) 
is the unperturbed energy of the bound state of the ground 
molecular potential. The photon of energy hfi is detuned 
from Eb\ by hSi, while the photon of energy /1/2 is detuned 
from Ef,2 by h&2- The decay rate of 61 is 71. Stark shifts 
of the levels due to trapping laser fields are neglected in this 
schematic. 



electric dipole-allowed transition due to the metastabil- 
ity of the 3 Pi state (r = 21.5/xs This increases 

the importance of the van der Waals interaction relative 
to the dipole-dipole term in determining the shape of 
the excited molecular potential, which makes the ground 
and excited potentials more similar than in alkali-metal 
atoms. Resulting Franck-Condon factors for transitions 
from free-atom to weakly bound excited molecular states 
are smaller, but overlap integrals between excited and 
ground molecular states are larger, which has implica- 
tions for the formation of ground state molecules, ft 
has been predicted that an optical Feshbach resonance 
induced by a laser tuned near an intercombination tran- 
sition can change the ground state scattering length 
with much lower inelastic loss 11811 than when using elec- 
tric dipole-allowed transitions |27l . [28| . Long coherence 
times are helpful for observing quantum optical effects, 
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FIG. 2: Atomic Sr energy levels involved in the two-photon 
PAS experiments (color online). Decay rates (s _1 ) and ex- 
citation wavelengths are given for selected transitions. Laser 
light used for the experiment is indicated by solid lines. Atoms 
decaying to the 3 P2 level may be repumped by 3 /im light. 



such as Autler-Townes splittings [29] of molecular levels 
[30l l3lL which is closely related to creation of an 
atom-molecule dark state @, HH and state-selective pro- 
duction of ultracold ground state molecules [3l], [U, [3|| . 



II. EXPERIMENTAL SETUP 

A. Laser Cooling and Trapping 

To perform two-photon spectroscopy, we start with 
laser-cooled atoms, and the initial cooling and trapping 
phases of the experiment are similar to previously pub- 
lished descriptions [Tl|, HH, Hi| . Atoms are trapped in a 
magneto-optical trap (MOT) operating on the 461 nm 
1 S'o- 1 Pi transition (Fig. ^ and cooled to about 2mK. 
There is a decay channel from the 1 Pi state, through the 
1 Z?2 state, to the metastable 3 p2 level with a branching 
ratio of 2 x 10~ 5 . To increase our sample number, we re- 
pump 3 P2 atoms by applying a 3 /im laser resonant with 
the 3 P2- 3 £>2 transition that returns these atoms to the 
ground state. The repumped sample of atoms contains 
about 3.5 x 10 s atoms. 

After this initial MOT stage, the 461 nm light is ex- 
tinguished and the atom sample is transferred with more 
than 50% efficiency to a second MOT operating on the 
1 5o- 3 Pi intercombination line @. The sample is cooled 
to 3/itK, producing densities of 10 12 cm~ 3 . The 689 nm 
light is provided by a master-slave diode laser system 
that is frequency-narrowed by servo-locking it to a high- 
finesse optical cavity with the Pound-Drever Hall method 
to produce a laser linewidth of ~ 50 kHz. Long-term sta- 
bility is maintained with a 1 Sq- 3 Pi saturated-absorption 
cell. 




FIG. 3: Photoassociation lasers (color online). The master 
laser that provides light for the intercombination-line MOT 
is frequency-stabilized via saturated absorption spectroscopy 
to the atomic transition, and it also provides the photoasso- 
ciation lasers. The one-photon PAS beam, with frequency /1, 
is generated directly from the master with an AOM. The two- 
photon PAS beam, with frequency /2, is formed by injection- 
locking a slave diode with a double-passed deflected beam 
from an AOM in a cat's eye configuration. 



B. Crossed Optical Dipole Trap 

To obtain high density and long sample lifetimes for 
improved two-photon PAS, atoms are transferred to an 
optical dipole trap (ODT) generated from a 21 W, 1064 
nm, linearly-polarized, single-transverse-mode fiber laser. 
The trap is in a crossed-beam configuration, derived from 
the first order deflection of an acousto-optic modulator 
(AOM). The beam is focused on the atoms with a min- 
imum e -2 intensity-radius of w=75 fim. It is then re- 
flected back through the chamber to intersect the first 
beam at 90 degrees and refocused to have the same waist 
at the atoms. Both beams lie in a plane that is inclined 
10.5 degrees from horizontal. The ODT trapping poten- 
tial is calculated from measured laser beam parameters 
and the polarizability of the 1 So state [37j • This allows 
us to determine the sample density profile from the tem- 
perature and number of trapped atoms. 

The maximum transfer efficiency observed from 
intercombination-line MOT to ODT for an optimized 
overlap time of 80 ms between the two is about 15%. This 
is limited in large part by inelastic collisions induced by 
689 nm light. Atoms are initially loaded with a single- 
beam ODT power of 5 W, which creates a trap depth of 
about U max /kB = 25 /iK. After the 689 nm light is extin- 
guished, the power is ramped to a final value between 2.5 
and 13 W in 20 ms, yielding equilibrium temperatures of 
between of 3 and 15 fiK. Up to 20 x 10 6 atoms are loaded 
to yield peak densities on the order of 10 14 cm -3 . 

The number of atoms and sample temperature are de- 
termined with time-of-flight absorption imaging using the 
1 S'o- 1 Pi transition. The lifetime of atoms in the ODT due 
to collisions with background atoms is about 2 s. 
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C. Photoassociation 

After the atoms have equilibrated in the ODT, the pho- 
toassociation lasers are applied (Fig.QJ. Laser fi is near 
resonance with a single-photon, free-bound transition to 
the red of the 1 So- 3 P\ atomic transition. For some stud- 
ies, this is the only laser applied. For two-photon PAS, 
laser ji is near resonance with a transition from the ex- 
cited molecular bound state to a ground molecular level. 

Using acousto-optic modulators, all photoassociation 
lasers are derived from the master laser that provides the 
intercombination-line MOT beams (Fig.[3|). Both fi and 
f% lasers are coupled into the same single-mode optical 
fiber with the same linear polarization. More than 85% 
of the intensity of both output beams is linearly polarized 
in the vertical direction, perpendicular to the ODT laser 
polarization. These beams are focused tow = 200 /im at 
the location of the atom sample, which is substantially 
larger than the atom cloud. The powers are monitored 
by a photodiode after the fiber. Alignment of the /i 
and /2 beams is facilitated by co-propagating the beams 
with 461 nm light aligned on the ODT atom sample by 
absorption imaging using an independent CCD camera. 
Depending on the specific measurement, one of the lasers 
is scanned and the number of atoms remaining after the 
photoassociative interaction time is recorded to obtain 
the loss spectrum. 



III. THEORETICAL DESCRIPTION OF 
PHOTOASSOCIATIVE LOSS 

Photoassociation is monitored by measuring the loss of 
ground-state atoms from the ODT. This loss is described 
with a local equation for the evolution of the atomic den- 
sity 



-2Kn 2 - Tn, 



(1) 



where the laser-frequency dependence of the collision 
event rate constant, K, determines the spectrum of the 
photoassociative loss. The observed PAS spectrum is rel- 
atively simple because the bosonic isotopes of strontium 
lack hyperfme structure. As shown in Fig. [I] ground state 
1 S atoms collide on a single 1 S+ potential. Four molec- 
ular potentials converge to the 1 S'o + 3 -Pi asymptote fl7j |. 
but only states of the U and l u potentials are optically 
excited from the 1 E+ potential |2(j. At the low temper- 
atures of atoms in the ODT, only s-wave collisions occur 
so only J = 1 intermediate levels and J = and 2 final 
states are populated. 

Photoassociative loss can be analyzed with the theory 
of Bohn and Julienne [H| , which yields 

K= hh : I \ S ( E ^f^h,...)?e- B ^ kBT dE g (2) 

where the partition function is Qt = O^jjr-^j f° r 
reduced mass y,. In spite of the low temperature, ther- 



mal averaging over the collision energy E g is necessary 
because of the narrow lincwidth of the transition. IS"! 2 
is the scattering probability for loss and its structure de- 
pends upon the loss process that is dominant. 

In two-photon spectroscopy of alkali-metal systems, 
the dominant photoassociative loss process is a collision 
on the open channel of two ground state atoms (g) with 
total energy E g leading to loss-producing decay from the 
excited state b\ with rate 71. (See Fig. [TJ) However, in 
the experiment reported here, b\ is metastable, and there 
is a concern that loss from the ground molecular state 62 
may also be important. The complete vanishing of the 
photoassociative loss when the lasers are on two-photon 
resonance from g — > 62 fSecs. fVl and IVI|) . however, im- 
plies that decay from 62 is negligible for our conditions 
and loss from b\ dominates. We can express the scatter- 
ing probability as 



I 



(3) 



(Si - S 2 ) 2 7 i7 s /(2tt) 2 



51 



5Sa 



} +(^) 2 (<5 1 -<5 2 f 



7i = diatomic, where 7 'atomic is the decay rate of the 
atomic 3 Pi level, and r ) s (E g ) is the stimulated width of 
bi due to laser-coupling to g, 
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where we represent g as the energy-normalized collid- 
ing state \E g ), leading to the the Franck-Condon fac- 
tor for the free-bound transition, | ( £> ! | £7^ ) | 2 , and V = 
( r \ 1/2 

( 2£^c ) f° r f rce- bound laser intensity I\ and molec- 
ular dipole matrix element d. Note that our 12 is the 
splitting of the Autler-Townes doublet (Sec. [Vj) , which 
differs from the Bohn- Julienne definition of the molecu- 
lar Rabi coupling [38|. 

The thermal energy is much greater than the zero- 
point energy for trap motion, T » hvtrap/kB, so con- 
finement effects are negligible [20]. We also neglect 
Doppler shift and photon recoil [16| . which is reason- 
able since T > Tr, where the recoil temperature for 
A = 689 nm photons is T R = h 2 /(k B \ 2 m) = 460 nK. 
We assume the decay products leave the trap, which is 
a good approximation for the intermediate levels we use, 
although it is not for the least bound U excited molecu- 
lar state (Mill. 

The energy integral for K (Eq. ^ is not analytic and 
must be evaluated numerically. The situation is further 
complicated by the ODT, which is not at the magic wave- 
length for one-photon photoassociation (914 nm [20|). 
The AC Stark shift of the weakly bound ground molec- 
ular level (62) is approximately equal to the shift of the 
incoming channel of 2-free atoms (g) [HI, H3]- In other 
words, the polarizability of the ground molecule is about 
twice that of a single atom. But the ground and excited 
molecular levels do not experience the same shift. For 
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spectroscopy, we can thus treat the ODT Stark shifts as 
a position-dependent shift of the intermediate state and 
dehne the laser detunings 



<*i = A - (E bl - E g )/h - X IoD T (r) 
h = h ~ (E bl - E b2 )/h - X IoDT{f), 



(5) 



where 1q dt {r) is the intensity profile of the optical dipole 
trap and x can be related to the differences in polariz- 
abilites for 1 So and 3 Pi atoms for 1.06 /jm laser light. 

This implies that |S| 2 and thus K are functions of po- 
sition, which must be addressed when Eq. [T]is integrated 
over the trap volume to calculate the time evolution of 
the number of trapped atoms 
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where Nq is the number at the beginning of the PAS 
interaction time. The one-body loss rate, T, is due to 
background collisions and off-resonant scattering from 
the PAS lasers. The effective volumes are defined by 



V q = I d 
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where U(r) is the trap potential, and 
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The kinetic energy integral is truncated by the local trap 
depth, U max — U(r). The spatial integrals in Eqs.[7]and 
[8] extend over the trap volume V in which U{r) < U m ax- 
Atom temperatures vary by no more than 25% during the 
interaction time, so assuming a constant sample temper- 
ature is reasonable. 

The spectrum is sensitive to many atomic and molec- 
ular parameters, and multiple types of spectra can be 
used to determine them. The ultimate goal is an accu- 
rate determination of E b2 because the molecular binding 
energy determines the s-wave scattering length and the 
underlying potential with high accuracy. 
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FIG. 4: Left: Atom number versus free-bound laser detuning 
from the one-photon 1 So- 3 Pi atomic transition (color online). 
Spectra shown here are for 6 W and 13 W ODT single-beam 
powers with sample temperatures of 6 fiK and 13 ^K, respec- 
tively. Right: Collision-event rate constant K e ff derived from 
the atom loss. The ODT at 1064 nm causes an AC Stark shift 
of the excited molecular state compared to the ground state, 
which shifts and broadens the line. The solid lines are fits 
using Eqs. [S]and[5] A peak shift of 480 kHz is measured for 
a single-beam power of 13 W. The dashed line marks the po- 
sition of our measured unperturbed resonance frequency at 
-222.25(15) MHz, which is in reasonable a gree ment with a 
previous measurement of -222.161(35) MHz [201 ]. 



The state iq is equal to the J = 1 rotational state of 
the third-least bound vibrational level of the 0„ poten- 
tial, with energy E bl = E 3Pl -hx 222.161(35) MHz [H- 
The Condon point for this excitation, where Vq u (R c ) — 
ViS (R c ) = E 0l , occurs at R c = 75 ao, which is very near 
the node in the ground state wave function [U [4l[ . 

Fitting data of atom number after a given interaction 
time, N(t), to Eq. O yields the collision-event rate con- 
stant K e ff. Figure [4] shows typical spectra for this tran- 
sition. 



IV. ONE-PHOTON PHOTOASSOCIATION 
A. One-Photon PAS Spectrum 

One photon PAS allows us to determine x, the relative 
light shift of states on the ground and excited molecular 
potentials (Eq. 0, and the stimulated width r ) s (E g ) (Eq. 

a. 

When 72 = and f2i2 = in Eq. [3] we recover the 
one-photon PAS scattering probability for loss through 
decay of b\ 
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B. Determining the Stimulated Width 

The area under the one-photon PAS line (Fig. [4]) can 
be related to molecular and experimental parameters 
through 
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Here, we have neglected truncation of the energy inte- 
gral, which is a small correction. The Wigner threshold 
law implies j s (E g ) cx I\^jE g . So for low laser inten- 
sity {j s (E g ) < 71), the expression A w (ls(E g )} is 
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FIG. 5: Area under one-photon PAS spectra versus free- 
bound laser intensity I\ (color online). The area can be re- 
lated to molecular and experimental parameters to determine 
the stimulated linewidth y 3 (E g ) of the PAS transition due to 
1%. For low free-bound laser intensities the area is indepen- 
dent of temperature and linearly dependent on 7i. 



independent of temperature and linearly dependent on 
intensity (Fig. [5]), where (• • •) refers to a thermal aver- 
age. The PAS saturation intensity for this transition 
(hat, pas l/T 1 / 2 ) defined as the intensity for which 
(ls(E g )) = 71, is 8 W/cm 2 for T = 13 ^K. Express- 
ing this in terms of an optical length for the transition 
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FIG. 6: Collision-event rate constant K e ff versus free-bound 
laser detuning from 1 Sq- 3 Pi atomic resonance for different 
bound-bound laser intensities (color online). These Autler- 
Townes doublets are measured with the bound-bound laser 
frequency fixed such that S2 « while scanning /1. The 
splitting of the spectra is given by the Rabi frequency f2i2 /2n 
and varies as y/Ti (Shown by the lines; spectra offset is pro- 
portional to y/h), where the bound-bound laser intensity I2 
is indicated in the legend. The asymmetry in the lineshapes 
arises from the bound-bound laser frequency being slightly off 
resonance from the bound-bound transition. The free-bound 
intensity Ii is constant for all four spectra at 0.05 W/cm 2 . 
The sample temperature is 8/1K. 



C. Modelling the Spectra and Determining the 
Relative AC Stark Shift 

Numerical integration of Eq. [8] to find K e f f using Eq. 
[9] for the scattering probability allows us to model the 
one-photon PAS spectra, and the relative light shift pa- 
rameter (Eq. [5]) can be varied to fit the data (Fig. [4]) . We 
find x = 160±30kHz/ (100 kW/cm 2 ), in good agreement 
with [33| , which yields a peak shift of 480 kHz for our 
deepest trap. The line shifts to the blue with more ODT 
laser intensity, showing the polarizability of 3 Pi atoms is 
less than the polarizability of atoms. 

At higher temperatures, the lineshapes in Fig. [4] pos- 
sess red tails, which result from the convolution of the 
Lorentzian spectrum with the Maxwcll-Boltzmann dis- 
tribution of collision energies [lj| [2(| ■ 



V. TWO-PHOTON AUTLER TOWNES 
SPECTRUM 

If PAS spectra are recorded in the same fashion as in 
Sec. lIVl except the bound-bound laser is added near reso- 
nance (t5 2 ~ 0) with a large intensity I 2 , the loss spectrum 
is modified due to the coupling between b\ and 62. This 



forms a A system and the line is split into an Autler- 
Townes doublet, with splitting given approximately by 
fii2/27r, where 



ni 2 =a|<6 1 |6 2 )| 7l J— . (11) 

y ^1 sat,atom 

The overlap integral is related to the Franck-Condon fac- 
tor, F, through F = |(6i|fo2)| 2 - The saturation inten- 
sity for the atomic 1 Sq — 3 Pi transition is I sa t,atom = 
7rfrc7 Qtom i C /(3A 3 ) = 3/iW/cm 2 . The rotational line 
strength factor, a, accounts for the change in dipole mo- 
ment from atom to molecule due to symmetry of wave 
function and projection on a rotating molecular axis 

Figure [6] shows several Autler-Townes spectra for vari- 
ous intensities I2 for a ground molecular state (b 2 ) equal 
to the J — rotational state of the least-bound ground 
vibrational level. This is the v = 62 level counting from 
the bottom of the well. We find that S1i 2 /2tt = 1 MHz for 
an intensity of .35 W/cm 2 , which yields F = 0.28 ± 0.06 
for a = The energy of the ground molecular 

state, Eb 2 , can be found from fits of data in Fig. but 
it is more accurately found by varying <5 2 with 8\ ~ 0, as 
we will discuss in Sec. IVII 
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FIG. 7: Collision-event rate constant K e ff versus frequency 
difference between free-bound and bound-bound lasers for 
spectroscopy of the J = 0, v — 62 level of the X X E^" poten- 
tial (color online). The free-bound laser frequency is fixed 
close to the one-photon PAS resonance and the intensity 
is h = 0.05 W/cm 2 . The bound-bound laser frequency is 
scanned, and its intensity is indicated in the legend. On 
two-photon resonance, PAS loss is suppressed due to quan- 
tum interference. The solid lines are fits using Eqs. [8] and [3] 
which yield E^/h — —136.7(2) MHz. The atom temperature 
is 8 /jK. 

The asymmetry in the line strengths in each doublet 
in Fig. [S] arises from the coupling-laser frequency be- 
ing slightly off resonance from the Stark-shifted bound- 
bound transition (82 7^ 0). But 82 is small and the scal- 
ing with intensity shows that the Autler-Townes splitting 
varies as v^2 as expected. 

These spectra show the potential of the system for 
quantum optics and ultracold molecule formation. The 
condition of no PAS loss when both lasers are on reso- 
nance has also been called a dark resonance Q, or an 
atom-molecule dark state [33| because the system has 
been put in an atom-molecule superposition state with 
vanishing excitation rate to b\. Such a state has also 
been proposed as a vehicle for creating large numbers of 
ground state molecules using STIRAP [42j]. The level of 
suppression near 8\ = shows the coherence of this su- 
perposition state in these experiments. It is noteworthy 
that the doublet is split by many linewidths even with 
moderate coupling-laser intensity because of the intrinsi- 
cally narrow spectrum of intercombination-line PAS. 



VI. TWO-PHOTON SUPPRESSION OF 
PHOTOASSOCIATION 

For determining the binding energy of molecular levels 
of the ground state potential, we hold the frequency of the 
free-bound laser fixed close to the one-photon resonance, 
Si ~ 0, and scan <5 2 - When #2 — £1 = 0, the system is in 
two-photon resonance from state g to 62 , and one-photon 
photoassociative loss is suppressed due to quantum inter- 
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FIG. 8: PAS suppression spectra for the J = 2, v = 62 level 
of the ground molecular potential, as described in Fig. [7] 
(color online). The sample temperature is 9 fj,K. The free- 
bound laser intensity is I\ — 0.04 W/cm 2 , and the bound- 
bound intensity is indicated in the legend. The solid line is a 
fit using Eqs. E] and El which yields E b2 /h = -66.6(2) MHz. 



Isotope 


V 


J 


Exp. 


Theory 


Diff. 


88 


62 





-136.7(2) 


-136.7 


0.0 


88 


62 


2 


-66.6(2) 


-66.5 


-0.1 



TABLE I: Observed ground molecular levels and experimental 
and theoretical level energies in frequency units (E^/h). 

ference. At this point, j\ — f% = (Eb2 — E g )/h, so the 
spectrum allows accurate determination of Eb 2 . An av- 
erage over Eg is necessary in order to properly account 
for thermal shifts of the resonance. 

Figure [7] shows a series of spectra taken at various 
bound-bound intensities for & 2 equal to the J = 0, v = 62 
state. Detuning of the free-bound laser frequency (/1) 
from the free-bound resonance, which depends on the 
ODT light shift (x) and collision energy (E g ), causes 
slight asymmetry in the lines and broadening, but this 
shape is reproduced with our model for K e ff (Eq. [5]) 
using parameters independently determined in previous 
sections. Since the initial and final states experience 
roughly equal light shifts due to the trapping laser, the 
ODT AC Stark shifts do not shift the resonance. No 
significant shift of the binding energy with laser power 
was observed, and we place an upper limit of 100 kHz for 
our highest intensity, I 2 — 0.55 W/cm 2 . We have also 
measured the binding energy of the J = 2, v = 62 state 
(Table []}, and a typical spectrum is shown in Fig. [8j 



VII. DETERMINATION OF THE SCATTERING 
LENGTH 

Binding energy measurements can be used to accu- 
rately determine the s-wave scattering length for the 
X 1 E+ ground electronic molecular potential and to ob- 
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tain information on the van der Waals coefficients, C„. 
For our analysis, the inner part (internuclear spacing 
R < 20.8 ao) of the potential is described with a recently 
published energy curve (43j derived from the Fourier- 
transform spectrum of Sr2 and additional information 
on the zero-energy ground-state scattering wave function 
from PAS [U [4l[ . A multipolar van der Waals expan- 
sion in C n /R n is used to represent the potential at longer 
range (R > 22.7 ao), and the gap between the two regions 
is bridged with a spline interpolation to insure a smooth 
connection. The wave functions are calculated using a 
full quantum calculation [2l| . 

As was the case in the analysis of (43|, we do not 
have enough information to independently determine all 
the van der Waals coefficients, and improve on the pre- 
cise relativistic many-body calculation of 44j , which gave 
C 6 = 3103(7) a.u., C 8 = 3.792(8) x 10 5 a.u. 



4.215 x 10 7 a.u. The last bound level (v 



and Cio = 
62, J = 0) 



is very extended, with its outer turning point at roughly 
R ~ 100 ao. At this point, the leading contribution to 
the dispersion energy, — Co/i? 6 — Cs/i? 8 — Cio/i? 10 , arises 
from the Cq term. In fact, using the values of C n above, 
the C% contribution is only 1.22% of that of Cq, while the 
Cio contribution accounts for only 0.02%. At a shorter 
distance R ~ 20 ao, similar to the maximum separation 
of states measured in [43| , these contributions are roughly 
30.5% and 8.5%, respectively. In (43[, the precise theo- 
retical values of [HJ for Cq and C§ were used to fit the 
value of Cio, since the energy levels (v — — 50) were 
more deeply bound and corresponded to shorter range 
than our levels (v = 62, J = and 2). We note that less 
precise values of Cq and C% from [45j were also considered 
in g|. 

If we use Cq as a fit parameter to match the binding 
energ y o f the J = 0, v = 62 state, assuming Cg and Cio 
from |44j. the best fit value is Cq = 3151(1) a.u. where 
the quoted uncertainty only reflects uncertainty in the 
measured binding energy (see Table Jj). The resulting 
88 Sr X : E+ s-wave scattering length is a — —2.0(2) ao, 
where the uncertainty also only reflects uncertainty in the 
measured binding energy. If we use the value of Cio = 
6.60 x 10 7 a.u. from |i| instead, we find Cq = 3116.0(5) 
a.u., in reasonable agreement with 44]. The resulting 
scattering length is a = —1.2(2) ao- If instead, we take 
Cq and Cg from [44[ , and fit Cio j as was done in [43| , we 
find Cio = 7.488 x 10 7 a.u. and a = -0.9(2) a (again, 
with the uncertainties reflecting the uncertainty in the 
measured binding energy). 

It is difficult to assess the uncertainties in a related 
to these coefficients. The most conservative assessment 
encompasses the full range of values quoted here; a = 
— 1.4(6) and Cq — 3130(20) a.u. The uncertainties in Cq 
and C% quoted in [44j are quite small, however, and no 
uncertainty is quoted for Cio- So that might give more 
credence to the results for Cq and C$ from [44j and the 
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2- phot. 
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1-phot. 


1-phot. 


Thermal- 
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transform 
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[46] 
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1000™ 


- 


430±|g 


84-84 
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- 
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1790(130) 
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< -1900 


- 
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- 


87-86 


162.5(5) 


165.5(5.5) 








87-84 


-56(1) 


-55(10) 








86-84 


31.9(3) 


33(3) 









resulting fit Cio = 7.488 x 10 a.u, which corresponds 
to the higher ends of the ranges of values for a in Table 
HT1 Mass-scaling can be used to determine the scatter- 



TABLE II: Published values of scattering lengths given in 
units of ao = 0.053 nm 



ing lengths for all stable-isotope collisional combinations 
from this information about the potential (Table [TTJ) . 

We note that the rotational energy takes the form 
H rot = B(R)J(J + 1) where B{R) = H 2 /{2fiR 2 ) is the 
rotational constant for separation R and reduced mass 
[i. The rotational constant for v = 62 calculated using 
this potential yields a binding energy of —66.5 MHz for 
the J = 2, v = 62 state (Table HJ, well within the mea- 
surement uncertainty. We also found that this quantity 
does not significantly constrain the C n coefficients. 

In the ultracold regime, the energy dependence of scat- 
tering lengths can often be neglected. However, this is 
not the case when there is a low-energy scattering reso- 
nance or when the zero-energy scattering length is very 
small. Figure [9] demonstrates that 88 Sr- 88 Sr and 86 Sr- 
86 Sr collision cross sections vary significantly with colli- 
sion energy, even at energies below 1 /iK. This may ex- 
plain the discrepancy between this work and previous 
studies of Sr collisional properties [46[ . 



VIII. CONCLUSION 

Using two-photon photoassociative spectroscopy, we 
have measured the binding energy of the J — and J = 2 
rotational levels of the v — 62 vibrational state of the 
X 1 S+ potential of 88 Sr2. This is the least-bound ground 
vibrational level. Combined with an accurate short range 
potential [43| and calculated van der Waals coefficients 
[44|, the measurement allows an accurate determination 
of a for 88 Sr- 88 Sr interactions. Through mass scaling, 
we determine the scattering lengths for all other isotopic 
combinations. These measurements serve as a stringent 
test of atomic structure calculations for alkaline-earth 
atoms and will provide valuable input for future experi- 
ments with ultracold strontium. 
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FIG. 9: Dependence of elastic-scattering cross sections on 
collision energy (E) in Kelvin for selected Sr isotopes (color 
online). The thick lines are cross sections including partial 
waves up to I = 4. Shape resonances are indicated. Thin 
lines indicate cross section contributions from 1 = only. For 
the plot, a potential is used that gives as8-88 = — 1.2ao at 
E — 0. The cross sections are given by the usual expressions: 

a = (8-K/k 2 ) J2Zo 2 ( 2l + D sin2 5l ( k ) k= >° 8?m2 for identical 

bosons and a = (4-ir/k 2 ) Y^ZoC 21 + 1 ) sin2 S i( k ) ^ 47 ™ 2 for 
distinguishable atoms. The phase shift, 5i=o(k), depends on 
k = y/2[iE/h and is related to the scattering length, a, and 
effective range, r e , at low k by k cot 5i-o(k) = — i + ^r e k 2 . 
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